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The jet shape modification due to partonic energy loss in the dense QCD matter is investigated 
by the help of the special transverse energy-energy correlator in the vicinity of maximum energy 
deposition of every event. In the accepted scenario with scattering of jet hard partons off comoving 
medium constituents this correlator is independent of the pseudorapidity position of a jet axis and 
becomes considerably broader (symmetrically over the pseudorapidity and the azimuthal angle) 
in comparison with pp-collisions. At scattering off "slow" medium constituents the broadening 
of correlation functions is dependent on the pseudorapidity position of a jet axis and increases 
noticeably in comparison with the previous scenario for jets with large enough pseudorapidities. 
These two considered scenarios result also in the different dependence of jet quenching on the 
pseudorapidity. 
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I. INTRODUCTION 

Experimental results from the CERN Super Proton 
Synchrotron (SPS) have revealed a rich set of new phe- 
nomena as the anomalous suppression of charmonium PD] 
and the broadening or mass shift of vector mesons |2( 
indicating the formation of dense matter 0, 0. How- 
ever, it has proven difficult to disentangle hadron con- 
tributions to the observed signals, and no clear consen- 
sus has emerged that a long lived quark-gluon plasma 
(QGP) has been observed in the fixed target heavy ion 
experiments at those facilities. Experiments at the BNL 
Relativistic Heavy Ion Collider (RHIC), which began op- 
eration in 2000, have initiated a new era in the study 
of QCD matter under extreme conditions. Nuclear col- 
lisions at the largest RHIC energy (^/snn = 200 GeV) 
not only produce matter at the highest energy density 
ever achieved but also provide a number of rare observ- 
ables that have not been accessible previously. Already 
in the first three years of RHIC operation the wealth 
of new data has been collected and analyzed such as a 
strong elliptic flow, the suppression of hadron spectra and 
azimuthal back-to-back two-particle correlations indicat- 
ing that a dense equilibrated system is generated in the 
most violent head-on collisions of heavy nuclei (see, for 
instance, a review [j| and references therein). 

The cm. energy for heavy ion collisions at the CERN 
Large Hadron Collider (LHC) = 5.5 TeV) will 

exceed that at RHIC by a factor of about 30. This pro- 
vides exciting opportunities for addressing unique physics 
issues in a new energy domain, where hard and semi-hard 
QCD multi-particle production can stand out against 
the underlying soft events. Besides the statistics is ex- 
pected to be high enough |fj for the systematic analysis 
of various aspects of QCD-physics (poorly accessible for 
studying at RHIC energy) in the medium with the ini- 
tial energy density considerably above the critical one 
for the deconfinement transition. In particular it will 
be possible to study the modification of jet shape and 



spectra in comparison with pp-collisions. Jets as one of 
the important probes of quark-gluon plasma are created 
at the very beginning of the collisions process (rf orm ~ 
1/pT < 0.01 fm/c) by the initial hard parton-parton scat- 
terings. These hard partons pass through the dense mat- 
ter formed at longer time scales (> 0.1 fm/c) and interact 
strongly with the medium constituents changing its origi- 
nal properties (the momentum direction, energy, particle 
distribution inside a jet) as a result of in-medium rescat- 
tcring. 

One of the encouraging methods for investigation of 
these effects seems to be the measurement of energy- 
energy correlations 0, H] ■ In the present paper we extend 
our approach 0, considering the energy-energy correla- 
tor in the restricted rapidity-angle space in the vicinity of 
maximum energy deposition of every event. This special 
correlator allows us to investigate, using the calorimct- 
ric information, the broadening of the jet shape due to 
partonic energy loss which is intensively discussed in the 
current literature IfJ. 



II. ENERGY CORRELATORS IN 
ELECTRON-POSITRON, HADRON AND 
NUCLEAR COLLISIONS 

To be clear at first we recall the well-known definition 
of the energy-energy correlation function £ which has 
been used by all for LEP experiments [nj at CERN and 
the SLD experiment [llj at SLAC to measure the strong 
coupling constant a s in e + e~-annihilation at the Z° res- 
onance with a high accuracy. £ is defined as a function of 
the angle x between two particles i and j in the following 
form: 

d ^(x) = o V V EiE i (l) 

dcos(x) Acos( X )iV event £^ E* ' {> 

where a is the total cross section for e + e~ — > hadrons, 
E is the total energy of the event, Ei and Ej are the 
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energies of the particles i and j. The sum runs over all 
pairs i, j with cos(xij) in a bin width Acos(x): 

cos(x) - Acos(x)/2 < cos(x 2 j) < cos(x) + Acos(x)/2. 

The limits A cos(x) — ► and N event — > oo have to be 
taken in Eq. QJ. 

This function can be calculated in perturbative QCD 
as a series in a s : 



1 dS(x) 
cr dcos(x) 
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(/3 ln(|)A( X ) + J B( X ))+O(a3), 



where [3o — (33 — 2n/)/3) "-/is the number of active 
flavors at the energy i£, Co is the Born cross section, \x 
is the renormalization scale. The first order term A(x) 
has been calculated by Basham et al. 12] from the well- 
known one gluon emission diagrams 7* , Z° — > qqg with 
the result 



A(x)= C F (l + u) 
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((2 - 6w 2 



xln(l + l/cj) + 6w-3), 
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where CV =4/3, oj = cot 2 (x)/2 and x is the angle be- 
tween any of the three partons initiating the three hadron 
jets. This allowed one to determine the strong coupling 
constant directly from a fit to the energy-energy correla- 
tions in e + e~-annihilation, since £ is proportional to a s 
in the first order. 

In hadronic and nuclear collisions jets are produced 
by hard scatterings of partons. In this case it is conve- 
nient to introduce the transverse energy-energy correla- 
tions which depend very weakly on the structure func- 
tions II3 and manifest directly the topology of events. 
Thus, for instance, in high-py two-jet events the corre- 
lation function is peaked at the azimuthal angle <p = 0° 
and ip — 180°, while for the isotropic background this 
function is independent of ip. 

Utilization of hard jet characteristics to investigate 
QGP in heavy ion collisions is complicated because of 
a huge multiplicity of "thermal" secondary particles in 
an event. Various estimations give from 1500 to 6000 
charged particles per rapidity unit in a central Pb+Pb 
collision at LHC energy, and jets can be really recon- 
structed against the background of energy flux beginning 
from some threshold jet transverse energy E^ ~ 50—100 
GeV 0, The transverse energy-energy correlation 

function is just sensitive to the jet quenching as well as to 
the global structure of energy flux (i.e. its anisotropy for 
non-central collisions) if we select events by an appropri- 
ate way 8] . The generalization of Eq. QJ for calorimetric 
measurements of the energy flow is straightforward: 
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FIG. 1: The region of r\ x <j> space used in order to calculate 
correlators in the near-side jet region. 



where 
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Ti 



is the total transverse energy in N calorimetric sectors 
covering the full azimuth, Et% is the transverse energy 
deposition in the sector i (i = l,...,N) in the consid- 
ered pseudorapidity region \r)\, n — [ip/Aip] (the integer 
part of the number ip/Aip), Aip ~ 0.1 is the typical az- 
imuthal size of a calorimetric sector. In the continuous 
limit Aip —> 0, N = [2ir/Aip] -> 00 Eq. Q reads 
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where E^ ls = 



dE T (4>) dE T 



,dE T {4>) dE T ((f> + ip) 



(5) 



is the distribution of 



transverse energy flow over azimuthal angle (f>. 

In our previous work Q it has been shown that at the 
special selection of events for the analysis (at least one 
high-px jet) and the procedure of event-by-event back- 
ground subtraction in high multiplicity heavy ion colli- 
sions, the transverse energy-energy correlator is sen- 
sitive to the partonic energy loss and angular spectrum 
of radiated gluons. The medium-modified energy-energy 
correlation function manifests significant strengthening 
in a wide interval of azimuthal angles around 7r/2, mod- 
erate broadening of the near-side jet region ip < 0.5 
and weak additional suppression of back-to-back corre- 
lations for ip ~ 7r. Without jet trigger this correlation 
function shows the global structure of transverse energy 
flux: the correlator is isotropic for central collisions and 
for non-central collisions it is sensitive to the azimuthal 
anisotropy of energy flow reproducing its Fourier har- 
monics but with the coefficients squared. 

Here we extend our approach 0, considering the 
energy-energy correlator in the restricted rapidity-angle 
space in the vicinity of maximum energy deposition of ev- 
ery event (Fig. 1). It is convenient to define this special 
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correlator in the following form: 
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where Eth is the transverse energy deposition in the 
pseudorapidity strip k (k = — M, — 1, 0, +1, +M) 
with the width A77 ~ 0.1 and the position of a jet axis in 
the strip k = 0, summarized over all azimuthal sectors i 
(i = -M, -1, 0, +1, +M) with the size A<p = Arj 
and the position of a jet axis in the angle sector i = 0; 
to = [ij/Af]} (the integer part of the number rj/Arf). Here 

E^ ls = ^2 Eti is the total transverse energy deposition 
i 

in the square of rj x <f> space, the center of which is deter- 
mined by the position of a jet axis. The similar definition 
can be written for the correlator as a function of the az- 
imuthal angle ( it is necessary to replace rj «-> <j> in Eq. 
© only) . The autocorrelation term £ Ej., is subtracted 

l 

in the denominator of the definition in order to normalize 
the integral (the sum) of correlator over -q to one in the 
fixed square around the jet with the number of pseudo- 
rapidity strips and angle sectors (2M + 1). 



III. A MODEL FOR ENERGY LOSS OF A JET 

We demonstrate the productivity and effectiveness of 
such correlators in the framework of our well worked-out 
model of jet passing through a medium. The model has 
been early applied to the calculation of various observ- 
ables sensitive to the partonic energy loss: the impact 
parameter dependence of the jet production |l5j , the 
mono/dijet rate enhancement [lfil J]J. the dijet rate de- 
pendence on the angular jet cone |18|. the elliptic coeffi- 
cient of the jet azimuthal anisotropy |19l l20l l2ll l22j , the 
anti-correlation between the softening jet fragmentation 
function and the suppression of the jet rate j2^]. The 
validity of approximations used has been verified by the 
comparison of RHIC data with the model calculations. 
The agreement is good j24j at the reasonable choice of 
main model parameters taking into account their depen- 
dence on the collision energy. At present the model has 
been constructed as the fast Monte Carlo event genera- 
tor PYQUEN (PYthia QUENched), and corresponding 
Fortran routine pyquen is available via Internet |25| . The 
routine is implemented as a modification of the standard 
PYTHIA jet events [2^. For details one can refer to 
these mentioned papers (mainly, 0,0, EE 13)- Here 
we note only the main steps essential for the present in- 
vestigation. 

The energy-energy correlator depends on not only the 
absolute value of partonic energy loss, but also on the 
angular spectrum of in-medium radiated gluons. Since 
coherent Landau-Pomeranchuk-Migdal radiation induces 
a strong dependence of the radiative energy loss of a jet 



on the angular cone size [18|,l22lL28, 29, 30], it will soften 
particle energy distributions inside the jet, increase the 
multiplicity of secondary particles, and to a lesser degree, 
affect the total jet energy. On the other hand, collisional 
energy loss turns out to be practically independent of the 
jet cone size and causes the loss of total jet energy, be- 
cause the bulk of "thermal" particles knocked out of the 
dense matter by elastic scatterings fly away in an almost 
transverse direction relative to the jet axis ^E]. Thus 
although the radiative energy loss of an energetic par- 
ton dominates over the collisional loss by up to an order 
of magnitude, the relative contribution of the collisional 
loss of a jet grows with increasing jet cone size due to 
the essentially different angular structure of loss for two 
mechanisms [l8|. Moreover, the total energy loss of a jet 
will be sensitive to the experimental capabilities to detect 
low-pT particles - products of soft gluon fragmentation: 
thresholds for a giving signal in calorimeters, influence of 
the strong magnetic field, etc. [3l|. 

Since the full treatment of the angular spectrum 
of emitted gluons is rather sophisticated and model- 
dependent [3 EE EE EE EH 5 we considered two sim- 
ple parameterizations of the distribution of in-medium 
radiated gluons over the emission angle 8. The "small- 
angular" radiation spectrum was parameterized in the 
form 
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where 6q ~ 5° is the typical ang le of the coherent gluon 
radiation estimated in Ref. 18]. The "broad-angular" 
spectrum has the form 



dN<> 1 

~dT a e 



(8) 



We believe that such a simplified treatment here is 
enough to demonstrate the sensitivity of the energy- 
energy correlator to the medium-induced partonic energy 
loss. 

PYTHIA_6.2 [2i| was used to generate the initial jet 
distributions in nucleon-nucleon sub-collisions at y/s — 
5.5 TeV. After that, event-by-event Monte Carlo simula- 
tion of rescattering and energy loss of jet partons in QGP 
was performed. The approach relies on accumulative en- 
ergy losses, when gluon radiation is associated with each 
scattering in the expanding medium together with in- 
cluding the interference effect by the modified radiation 
spectrum as a function of decreasing temperature. Such 
a numerical simulation of the free path of a hard jet in 
QGP allows any kinematical characteristic distributions 
of jets in the final state to be obtained. Besides, the 
different scenarios of medium evolution can be consid- 
ered. In each ith scattering a fast parton loses energy 
collisionally and radiatively, Ae^ = tj/(2mo) +ct>», where 
the transfer momentum squared tj is simulated accord- 
ing to the differential cross section for elastic scattering of 
a parton with energy E off the "thermal" partons with 
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energy (or effective mass) m ~ 3T <C E at tempera- 
ture T, and u>i is simulated according to the energy spec- 
trum of coherent medium-induced gluon radiation in the 
Baier-Dokshitzer-Mueller-Schiff formalism (27). Finally 
we suppose that in every event the energy of an initial 
parton decreases by the value AE = J^i Ae,. 

The medium was treated as a boost-invariant longi- 
tudinally expanding quark-gluon fluid, and partons as 
b eing produced on a hyper-surface of equal proper times 
r |32j . For certainty we used the initial conditions for the 
gluon-dominated plasma formation expected for central 
Pb+Pb collisions at LHC H|: r ~ 0.1 fm/c, T ~ 1 
GeV. For non-central collisions we suggest the propor- 
tionality of the initial energy density to the ratio of the 
nuclear overlap function and the effective transverse area 
of nuclear overlapping 15]. 

In the frame of this model and using above QGP 
parameters we evaluate the mean energy loss of quark 
of Et = 100 GeV in minimum-bias Pb+Pb collisions, 
< AEj, >~ 25 GeV (accordingly, the gluon energy loss 
is < AiSf, >= 9/4 < AE^ >). In our simulation the 
value of energy loss in fact is regulated mainly by the 
initial temperature Tq. 



IV. NUMERICAL RESULTS AND DISCUSSION 

To be specific, we consider the geometry of the Com- 
pact Muon Solenoid (CMS) detector [3l| at LHC. The 
central (barrel) part of the CMS calorimetric system cov- 
ers the pseudorapidity region \rj\ < 1.5, the segmenta- 
tion of electromagnetic and hadron calorimeters being 
Ar/x A(f) = 0.0174x0.0174and ArjxAcf) = 0.0872x0.0872 
respectively • The endcap part of the CMS calorime- 
ter 1.5 < \rj\ < 3 has more tangled structure: the tower 
sizes of hadron calorimeter over r) and 4> are 0.087 at 
1 77 1 = 1.5 and increase gradually up to 0.345 at |ry| = 3. 

Then the energy-energy correlation function JSJ is cal- 
culated for the events containing at least one jet with 
E*t > E T (threshold) = 100 GeV in the considered kinc- 
matical region. The final jet energy is defined here as 
the total transverse energy of final particles collected 
around the direction of a leading particle inside the cone 
R = \J Arf 2 + A(f> 2 = 0.7, where 77 and <fi are the pseudo- 
rapidity and the azimuthal angle respectively. Note that 
the estimated event rate in the CMS acceptance, ~ 10 7 
jets with Et> 100 GeV in a one month LHC run with 
lead beams ja, |31| , will be large enough to carefully study 
the energy-energy correlations over the whole azimuthal 
range. 

We restrict ourself to central collisions in which the 
effect of jet quenching is maximum and take into con- 
sideration the energy deposition resulting from jets only. 
The fact is that in the CMS heavy ion physics program 
the modified sliding window-type jet finding algorithm 
has been developed to search for "jet-like" clusters above 
the average energy and to subtract the background from 
the underlying event [(J 0, This developed algo- 



rithm allows one to reconstruct the jet energy and its 
space position with a high enough accuracy in a high 
multiplicity environment beginning from some threshold 
jet transverse energy Ej. ~ 50 — 100 GeV. Thus, for 
instance, the reconstruction efficiency (the ratio of the 
number of reconstructed jets to the total number of gen- 
erated jets) for jets with energy E^ > 100 GeV is close 
to 100%, the fraction of false jets being not greater than 
1% of the total number of reconstructed jets. The mean 
reconstructed jet energy E^fceco) is a linear function of 
the mean generated jet energy By'(MC), the shape of 
this dependence being identical in the presence and in 
the absence of the background from Pb+Pb events. This 
means that pp-collisions can be used for an adequate com- 
parison with Pb+Pb collisions. The attainable precision 
of the space resolution of jets is rather high: it is less 
than the size of a hadron calorimeter tower 0, l3lL l34| . 

This gives grounds to hope that in the region of max- 
imum energy deposition the contributions from jets and 
"thermal" background can be effectively separated not 
only on the simulation level, when they are known a pri- 
ori. Therefore in order to reveal the new possible issues 
of the investigation of energy correlations we neglect this 
possible difference between the initially generated and 
reconstructed jet energy depositions, which can be sig- 
nificant on the outlying jet area only. Though the analy- 
sis of the quality of extracting the jet-like energy-energy 
correlator against background fluctuations based on the 
detailed simulation of detector responses (together with 
the optimization of the jet detection threshold and the 
calibration coefficients used at the subtraction and re- 
construction procedures) should be performed for each 
specific experiment. In particular, in CMS we have to 
take into consideration that most of the charged low-pr 
particles will be cleared out of the central calorimeters by 
the strong magnetic field (decreasing a really measurable 
energy flux and absolute value of its fluctuations in such 
a way) . 

The results of the numerical simulations are presented 
in Fig. 2 (for the cases without and with partonic energy 
loss) for the events containing at least one hard jet with 
E^ > £+ (threshold) = 100 GeV. The two parameteriza- 
tions of the distribution on the gluon emission angles Q 
and © were used and At? = Acj) = Atp = 0.087, M = 8. 
One can see that the energy loss results in the noticeable 
broadening of jet shape which can be characterized by 
the variation of the mean pseudorapidity squared (the 
width of distribution) . These mean values of pseudora- 
pidity squared are defined in the following form: 

<rf> = J^^^/Jd,-^ (9) 

and are displayed in Table I. They are mainly determined 
by the first bins of histograms closed to a jet axis where 
the separation of jet and background contributions to the 
total energy deposition is expected to be effective. The 
broadening effect is considerable for the "broad-angular" 
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FIG. 2: The transverse energy-energy correlator dTir/dn in 
the near-side jet region as a function of pseudorapidity r\ 
without (solid histogram) and with medium-induced partonic 
energy loss for the "small-angular" J7J (dashed histogram) 
and the "broad-angular" JSj (dotted histogram) parameteri- 
zations of emitted gluon spectrum in central Pb+Pb collisions 
for a jet axis position from the central pseudorapidity region 
|r? jct |<0.8. 



FIG. 3: The transverse energy-energy correlator dY^r/dip in 
the near-side jet region as a function of azimuthal angle ip 
without (solid histogram) and with medium-induced partonic 
energy loss for the "small-angular" J7J (dashed histogram) 
and the "broad-angular" JSJ (dotted histogram) parameteri- 
zations of emitted gluon spectrum in central Pb+Pb collisions 
for a jet axis position from the central pseudorapidity region 
|77 jct |<0.8. 



radiation, comes to ~ 45% and probably can be observ- 
able in spite of the ambiguities discussed above. 

Figure 3 shows the behavior of correlator ijfJJl in the 
near-side jet region as a function of azimuthal angle ip 
(the notations and the kinematical restrictions are the 
same as in Fig. 2). As in the previous case the broad- 
ening is more pronounced for the "broad- angular" pa- 
rameterization of emitted gluon spectrum. Besides the 
equality (with a good accuracy) of mean values of pseu- 
dorapidity and azimuthal angle squared results from the 
rapidity-angle symmetry of jet shape for jets from the 
central pseudorapidity region. For the illustration the 
difference between correlators 



dip 



v=<p 



as a function of r\ = ip is presented in Fig. 4. This 
difference is small, nonregular and comes to the level of 
statistical fluctuations. 

We are interested in the nonsymmetrical behavior of 
correlators as functions of pseudorapidity and azimuthal 
angle because there is some discussion in the current liter- 
ature |ij about the possible nonsymmetrical modification 
of jet shape due to the collective motion of medium con- 
stituents, on which the rescattering of jet partons takes 
place. The correlators under consideration could show 



the different behavior over angle and pseudorapidity for 
a jet axis position from the nonsymmetrical pseudorapid- 
ity region. However our simulation gives practically the 
former results for correlators calculated for jets from the 
noncentral pseudorapidity region 0.8 < n ict < 2.4. The 
corresponding values of pseudorapidity and angle squared 
are listed in the lower part of Table I. The difference be- 
tween correlators as a function of r] = ip is again small, 
nonregular and comes to the level of statistical fluctua- 
tions as in Fig. 4. 

Thus the behavior of transverse energy-energy corre- 
lators in the near-side jet region is independent of the 
pseudorapidity position of a jet axis. The fact is that in 
our model the rescattering of a hard parton takes place 
off a comoving medium constituent (i.e. moving with the 
same longitudinal rapidity as a hard parton). In such 
approach gluons are emitted at the polar angle 9 relative 
to the transverse momentum px of a radiating parton 
isotropically over the azimuthal angle in the rest frame 
of a medium constituent that, as a matter of fact, pro- 
vides for the longitudinal boost invariance of correlators. 
However this scenario with the same longitudinal rapid- 
ity of a jet parton and a medium constituent (mostly used 
and closely related with the Bjorken scaling solution) is 
not unique and one can suppose that a jet parton and 
a medium constituent move in the longitudinal direction 
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TABLE I: The mean values of pseudorapidity and angle squared as a measure of the width of transverse energy-energy correlation 
functions. 
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FIG. 4 
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The difference between correlators 



dr] 



in the near-side jet region as a function of r\ = ip 

without (solid histogram) and with medium-induced partonic 
energy loss for the "small-angular" J7J (dashed histogram) 
and the "broad-angular" J5J (dotted histogram) parameteri- 
zations of emitted gluon spectrum in central Pb+Pb collisions 
for a jet axis position from the central pseudorapidity region 
|77 jot |<0.8. 



in a different way. 

On the generator level it is not difficult to realize the 
scenario when the rescattering takes place off a so-called 
"slow" medium constituent (i.e. moving with the zero 
or considerably less longitudinal rapidity in comparison 
with the rapidity of a jet parton) and gluons are emitted 
at the polar angle 6 relative to the total momentum p of 
a parent parton. From Table I one can see that for jets 
from the central pseudorapidity region \r]' 3et \ < 0.8 and 
the "small- angular" parameterization J7J this modifica- 
tion does not lead to any variation of correlation func- 



tions because of the small distinction between the total 
and transverse momenta in this case. For the "broad- 
angular" radiation (JSJ the broadening of jet shape is 
independent of the rescattering scenario in the central 
pseudorapidity region at the additional restriction on the 
difference between the directions of an initial hard par- 
ton and a final jet axis (due to the large "tail" in the 
distribution over this difference). 

For jets from the forward pseudorapidity region 0.8 < 
jyjct < 2.4 (where the distinction between the total and 
transverse momenta is significant) the broadening of jet 
shape increases considerably that is supported by the cal- 
culated mean values of pseudorapidity and angle squared. 
This larger broadening results from that the energy loss 
is determined by not the transverse momentum (as this 
was early) but the total one which is considerably larger 
than the transverse momentum for jets from the forward 
pseudorapidity region. The experimental observation of 
such dependence of the width of correlation functions on 
the pseudorapidity position of a jet axis will indicate the 
scenario with the gluon radiation around the total mo- 
mentum of a parent parton. Besides this observation does 
not demand the comparison with pp-collisions, in which 
such dependence is absent. 

As regards the possible 77 x p asymmetry, then it is in- 
visible as before on the level of the mean values, although 
we might catch some regular (but again practically on the 
level of statistical fluctuations) difference in the behav- 
ior of correlators over rj and tp. This can mean that the 
transverse energy-energy correlations (unlike the average 
jet energy and jet multiplicity 0) is weakly sensitive to 
longitudinal flow effects due to the longitudinal Lorentz 
invariance of transverse energy itself or our simple im- 
itation of flow effects on the generator level can be, in 
principle, not enough to their study in respect of the 
possible 77 x ip asymmetry. 

The scenarios suggested above result also in the differ- 
ent dependence of jet quenching on the pseudorapidity. 
In the accepted scenario (I) with the rescattering off co- 
moving medium constituents jet quenching is expected 
to be independent of the jet pseudorapidity under the 
condition that the matter density is also ^-independent 
and Figure 5 demonstrates this independence. We define 
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FIG. 5: The jet quenching factor R as a function of pseu- 
dorapidity r\ without (solid line) and with medium-induced 
partonic energy loss for the "small-angular" Q (dashed his- 
togram) and the "broad-angular" © (solid histogram) pa- 
rameterizations of emitted gluon spectrum in central Pb+Pb 
collisions in scenario I. 

a jet quenching factor R in a conventional way, 

N(E iet > 100 GeV)\ with losses /1 n\ 

N(E ict > 100 GeV)| 

without losses 

as the ratio of the jet number with transverse energy 
E^ > 100 GeV with energy loss to the corresponding 
number of jets without loss (i.e. in pp-collisions nor- 
malized on the number of binary nucleon-nucleon sub- 
collisions). In the second scenario (II) a quenching fac- 
tor is strongly (Fig. 6) dependent on the pseudorapidity 
position of a jet axis. For the "broad-angular" radia- 
tion this factor R has more tangled behavior with the 
shifted (not in the zero) maximum over rj. This can be 
explained by the large "tail" in the distribution over the 
difference (especially in 77 X ip space !) between the di- 
rections of an initial hard parton and a final jet axis and 
by our simplified jet definition. Here the final jet energy 
is defined simply as the total transverse energy of final 
particles collected around the direction of a leading par- 
ticle inside the cone R = ^Arj 2 + Acj) 2 = 0.7 without 
correction of the jet axis usually used in the modified 
sliding window-type jet finding algorithm. The structure 
with the pseudorapidity-shifted maximum disappears if 
the restriction on the difference between the directions of 
an initial hard parton and a final jet axis is established as 
one can see in Fig. 7 or if the modified sliding window- 
type jet finding algorithm [(| I2J 0| is used. But for 
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FIG. 6: The jet quenching factor R as a function of pseu- 
dorapidity r\ without (solid line) and with medium-induced 
partonic energy loss for the "small-angular" Q (dashed his- 
togram) and the "broad-angular" © (solid histogram) pa- 
ramctcrizations of emitted gluon spectrum in central Pb+Pb 
collisions in scenario II. 



all that the additional restriction and the use of other 
algorithm have no influence on the rapidity-behavior of 
quenching factor in the scenario I. Thus the quenching 
factor together with the jet shape broadening is inde- 
pendent of the pseudorapidity position of a jet axis in 
the scenario I while in the scenario II the dependence 
of these observables on the pseudorapidity is noticeable 
and probably can be experimentally detected in spite of 
ambiguities discussed above. 



V. SUMMARY 

The special correlator in the vicinity of maximum en- 
ergy deposition of every event allowed us to investigate 
the jet shape modification due to partonic energy loss us- 
ing the calorimetric information. In the accepted scenario 
with scattering of jet hard partons off comoving medium 
constituents this correlator is independent of the pseudo- 
rapidity position of a jet axis and becomes considerably 
broader (symmetrically over the pseudorapidity and the 
azimuthal angle) in comparison with pp-collisions. At 
scattering off "slow" medium constituents the broaden- 
ing of correlation functions is dependent on the pseudo- 
rapidity position of a jet axis and increases noticeably in 
comparison with the previous scenario for jets with large 
enough pseudorapidities. These two considered scenarios 
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result also in the different dependence of jet quenching 
on the pseudorapidity. We believe that such a transverse 
energy-energy correlation analysis may be useful at LHC 
data processing. 
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